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Abstract 

The most general reaction-diffusion model on a Cayley tree with nearest- 
neighbor interactions is introduced, which can be solved exactly through 
the empty-interval method. The stationary solutions of such models, as 
well as their dynamics, are discussed. Concerning the dynamics, the spec- 
trum of the evolution Hamiltonian is found and shown to be discrete, 
hence there is a finite relaxation time in the evolution of the system to- 
wards its stationary state. 
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1 Introduction 



Reaction-diffusion systems have been studied using various methods, including 
analytical techniques, approximation methods, and simulation. Approximation 
methods are generally different in different dimensions, as for example the mean 
field techniques, working good for high dimensions, generally do not give cor- 
rect results for low-dimensional systems. A large fraction of analytical studies 
belong to low-dimensional (specially one-dimensional) systems, as solving low- 
dimensional systems should in principle be easier. [1-11]. 

The Cayley tree is a tree (a lattice having no loops) where every site is 
connected to £ nearest neighbor sites. This no-loops property may allow ex- 
act solvability for some models, for general coordination number £. Reaction 
diffusion models on the Cayley tree have been studied in, for example [12-17]. 
In [12,13,16] diffusion-limited aggregations, and in [14] two-particle annihilation 
reactions for immobile reactants have been studied. There are also some exact 
results for deposition processes on the Bethe lattice [17]. 

The empty interval method (EIM) has been used to analyze the one dimen- 
sional dynamics of diffusion-limited coalescence [18-21]. Using this method, 
the probability that n consecutive sites are empty has been calculated. This 
method has been used to study a reaction-diffusion process with three-site in- 
teractions [22] . EIM has been also generalized to study the kinetics of the g-state 
one-dimensional Potts model in the zero-temperature limit [23]. In [18-21], one- 
dimensional diffusion-limited processes have been studied using EIM. There, 
some of the reaction rates have been taken infinite, and the models have been 
worked out on continuum. For the cases of finite reaction-rates, some approxi- 
mate solutions have been obtained. 

In [24,25], all the one dimensional reaction-diffusion models with nearest 
neighbor interactions which can be exactly solved by EIM have been found and 
studied. Conditions have been obtained for the systems with finite reaction 
rates to be solvable via EIM, and then the equations of EIM have been solved. 
In [25], general conditions were obtained for a single-species reaction-diffusion 
system with nearest neighbor interactions, to be solvable through EIM. Here 
solvability means that evolution equation for E n (the probability that n con- 
secutive sites be empty) is closed. It turned out there, that certain relations 
between the reaction rates are needed, so that the system is solvable via EIM. 
The evolution equation of E n is a recursive equation in terms of n, and is linear. 
It was shown that if certain reactions are absent, namely reactions that produce 
particles in two adjacent empty sites, the coefficients of the empty intervals in 
the evolution equation of the empty intervals are n-independent, so that the 
evolution equation can be easily solved. The criteria for solvability, and the so- 
lution of the empty-interval equation were generalized to cases of multi-species 
systems and multi-site interactions in [26-28]. 

In this article the most general single-species reaction-diffusion model with 
nearest-neighbor interactions on a Cayley tree is investigated, which can be 
solved exactly through the empty interval method. The scheme of the paper is 
as follows. In section 2, the most general reaction-diffusion model with nearest- 
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Figure 1: The Cayley tree with £ = 3 



neighbor interactions on a Cayley tree is studied, which can be solved exactly 
through EIM. The evolution equation of E n is also obtained. In section 3 the 
stationary solution of such models, as well as their dynamics are discussed. 
Finally, section 4 is devoted to concluding remarks. 

2 Models solvable through the empty interval 
method on a Cayley tree 

The Cayley tree is a tree (a lattice without loops) where each site is connected 
to £ sites (fig. 1). Two sites are called neighbors iff they are connected through 
a link. Consider a system of particles on a Cayley tree. Each site is either 
empty or occupied by one particle. The interaction (of particles and vacancies) 
is nearest neighbor. The probability that a connected collection of n sites be 
empty is denoted by E n . It is assumed that this quantity does not depend on 
the choice of the collection. An example is a tree where the probability that a 
site is occupied is p and is independent of the states of other sites. Then 

E n = (l-p) n . (1) 

The following graphical representations help express various relations in a more 
compact form. An empty (occupied) site is denoted by o (•). A connected 
collection of n empty sites is denoted by On ■ 

There is no loop in a Cayley tree, so each site can only be connected to a 
single existing cluster site, by a single link. For £ > 3 (the case we are interested 
in here) the closedness of the evolution equation for E n requires that the rate 
of creating an empty site be zero. The reason is that if it is not the case, then 
an empty n-cluster can be created from two disjoint empty clusters joined by a 
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Figure 2: An empty cluster with the links at the boundary, on a Cay ley tree 
with £ = 3 



single occupied site [29] . This shows that if the evolution of the empty clusters 
is to be closed, then the only possible reactions are the following, with the rates 
indicated. 

•o — y • •, r\ 

OO — > O •, T2 

oo-»««, r 3 . (2) 

(There is no distinction between left and right, of course.) This means that the 
reactants are immobile, and the coagulation and diffusion rates are zero. 
Using these, one arrives at the following time evolution for E n : 

dE 

-^ = -R n ri P(-O n ) - R n (r 2 + r 3 ) P(o-O n ) 



-(n-l)(2r 2 + r 3 )P(OJ, (3) 

where i? n is the number of sites adjacent to a collection of n connected sites. A 
simple induction shows that 

R n = n (£ - 2) + 2. (4) 

One has 

P(*On) + P(°-On) = P (On), (5) 

from which 

P(»O n ) = E n - E n+1 . (6) 
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Using this, one arrives at 



dE 

— f - Rn [-n (£„ - £„+i) - (r 2 + r 3 ) £ n+1 ] - (n - 1) (2r 2 + r 3 ) £„. (7) 
Throughout the paper, it is assumed that n, ra, and r 3 are all nonzero. 

3 The solution 

The stationary solution of the system (_E S , for which the time derivative van- 
ishes), satisfies 

R n [-n (££ - £* +1 ) - (ra + r 3 ) E s n+1 ] - (n- 1) (2r 2 + r 3 ) = 0. (8) 

As En s are nonnegative and nonincreasing in n, it is easy to see that the only 
solution to © is 

K = 0- (9) 

This means that in the stationary configuration, all of the sites are occupied, 
which is not a surprise since in all reactions particles are created. 

Regarding dynamics, one question is to obtain the spectrum of the evolu- 
tion Hamiltonian. This is equivalent to finding solutions with exponential time 
dependence: 

= .E£exp(£t). (10) 

Putting this in (J7J), one arrives at 

- [Rn n + (n - 1) (2 r 2 + r 3 ) + £} E £ n + R n {n - r 2 - r 3 ) Ef l+l = 0. (11) 



From this, 
where 

It is seen that 



E n+l = CnE n , (12) 

, _ R n r 1 + {n-l)(2r 2 +r 3 )+£ 
R n (ri - ra - r 3 ) 

(g-2)n+2r 2 + r 3 

lim C« = 77 — ^-7 7- (14) 

n-»oo (f - 2) (n - ra - r 3 ) 

The right-hand side is either negative or greater than one. So if all E^s are 
nonzero, then E^'s either are not all nonnegative or blow up for large n's. Such 
E^s are not acceptable as probabilities. To see the reason, consider £\ (the 
largest £). for large times, only E n J s corresponding to this eigenvalue survive. 
But these should be nonincreasing with respect to n, and nonnegative, which is 
not the case. So E^ 1 's must be identically zero for n larger than a certain integer 
(say ?ii). A similar reasoning can then be made for £ 2 (the next largest value 
of £), and the values of E^ 2 for n > m, to show that there should be another 
integer n 2 so that E^ 2 vanishes for n > n 2 - This argument can be continued 
to show that for all £'s, there must be an integer so that E £ 's are identically 
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zero for n larger than that integer. This shows that ( n must be zero for some 
positive n, which gives the allowed values of £ : 

£ k = -Zn-{k-i)0, k>i, (15) 

where 

/3:=(£-2)r 1 +2r 2 + r 3 . (16) 

This spectrum is discrete, and there is a gap between the largest eigenvalue and 
zero, which means that the system evolves towards its stationary configuration 
with a relaxation time. This relaxation time is 

r=^~- (17) 

One can also find E%'s. Denoting E^ k by E%, and using ([12]) and ([IB"]), one 
arrives at 

E n = ^ > (18) 



(19) 



T\n+—j (k-n)\ 

where 

(£-2) (r 2 +r 3 - ri ) 
(f - 2) n + 2 r 2 + r 3 ' 
The general solution to © is then 



^n(t) = J2 Ck E n e M£k t), (20) 

fc=l 

where c^'s are to be determined from the initial condition. 
A special solution to (J7J) is of the form 

£ n (t)=S 1 (t)[6(t)]"- 1 . (21) 

Putting this in ([7]), one arrives at 

'. = -pb-{3ab 2 , 
at 

— + eh"") *■ (22) 

These are readily solved and one obtains 
&(0)exp(-/3t) 



&(*) = 



l + ab(0) [1 -exp(-0t)Y 



=^'»'-i'<-<n n { 1 + am[1 \^ { _ m y~' ■ (23) 



i 
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Using these, one obtains 

E n {t) = E n {0) exp[-£n t-{n-l)(3t] 



l + a6(0) [l-exp(-/3t)] 



(24) 



It is seen that for large times, all £? n 's tend to zero. In fact they decay like 

E n (t)~ exp[-Znt-(n-i)pt]. (25) 

One notes that in fact E n (t) decays like exp(— £ n t), and this is expected, as _E„ 
is zero for k < n. 

A special case where the ansatz ([21]) works is the case of initially uncorrelated- 
sites, so that each site is occupied with probability p regardless of other sites. 
One has then 

E n (0) = (l-p) n , (26) 



so that 



E 1 (p)=l-p, 
6(0) =1 - p. 



(27) 



The special case £ = 2 can be treated directly or as a limiting case of the 
general problem. The results corresponding to (fTS"]) and (fT5|) would be 



-2n-(fc-l)(2r 2 + r 3 ), £ 



and 



1 



(k-n)\ 



2 (r 2 + r 3 - ri) 
2 r 2 + r 3 



k—n 



£ = 2. 



Finally, the solutions corresponding to the ansatz (f2"Tj) would be 
6(f) =6(0) exp[-(2r 2 + r 3 )i], £ = 2, 



(28) 
(29) 

(30) 



and 



Ei(t) = Si(0) exp 



2 (n - r 2 - r 3 ) 



2r 2 +r 3 
xexp(-2rit), £ = 2, 



6(0) l-exp[-(2r 2 +r 3 )t] 



(31) 



so that 



= E n (0) exp 



2 (n - r 2 - r 3 ) 



2 r 2 + r 3 

x exp{-[2n + (n - 1) (2r 2 +r 3 )]i}, 



6(0) 1 -exp[-(2r 2 +r 3 )i] 
£ = 2, 



(32) 



G 



4 Concluding remarks 



The most general single-species exclusion model on a Cayley tree was consid- 
ered, for which the evolution of the empty-intervals is closed. It was shown 
that in the stationary configuration of such models all sites are occupied. The 
dynamics of such systems were also studied and it was shown that the spectrum 
of the evolution Hamiltonian is discrete. The time evolution of the initially 
uncorrelated system was also obtained. Among the questions remaining, one 
can mention the problem of Cayley trees with boundaries, with injection and 
extraction at the boundaries. 
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